Diffuse Neutron Scattering Study of a Disordered Complex Perovskite 

Pb(Zn 1/3 Nb 2/3 )03 Crystal 



o 
o 

(N 

Cu- 
te 

m 



D. La-Orauttapong, 1 J. Toulouse, 1 J. L. Robertson, 2 and Z.-G. Ye 3 

1 Physics Department, Lehigh University, Bethlehem, PA 18015-3182 
2 Oak Ridge National Laboratory, Solid State Division, Oak Ridge, TN 37831-6393 
3 Department of Chemistry, Simon Fraser University, Burnaby, BC, V5A 1S6, Canada 

Diffuse scattering around the (110) reciprocal lattice point has been investigated by elastic neutron 
scattering in the paraelectric and the relaxor phases of the disordered complex perovskite crystal- 
Pb(Zn 1; /3Nb2/3)03 (PZN). The appearance of a diffuse intensity peak indicates the formation of 
polar nanoregions at temperature T* , approximately 40K above T C =413K. The analysis of this 
diffuse scattering indicates that these regions are in the shape of ellipsoids, more extended in the 
< 111 > direction than in the < 001 > direction. The quantitative analysis provides an estimate 
of the correlation length, £, or size of the regions and shows that f/m) ~ 1.2f/ oi), consistent with 
the primary or dominant displacement of Pb leading to the low temperature rhombohedral phase. 
Both the appearance of the polar regions at T* and the structural transition at T c are marked by 
kinks in the £(m) curve but not in the £(ooi> one i a l so indicating that the primary changes take 
place in a < 111 > direction at both temperatures. 
PACS Numbers: 77.84.Dy,61.12.Ld,77.80.-e 
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Many of the relaxor ferroelectrics known today are 
lead-based compounds with the perovskite structure. In 
addition to the characteristic frequency dispersion of 
their dielectric constant, several of them exhibit remark- 
able piezoelectric or electrostrictive properties that are 
finding important applications, e.g. as transducers and 
actuators. The Pb 2+ -containing relaxor perovskites such 
as Pb(R 1/3 Nb 2 /3)0 3 (R=Mg 2 +, Zn 2 +) have a common 
ABO3 cubic perovskite structure in which the B-site can 
be occupied by |R 2+ and |Nb 5+ . Because of the dif- 
ferent atomic radii and valences of the B-site cation, 
PMN and PZN exhibit short-range chemical ordering 
Over the years, these relaxors and their proper- 
ties have been described in a variety of ways, most often 
in terms of the formation of polar micro- or nanoregions 
HQ. However, these models are primarily based on in- 
direct experimental evidence for such regions and more 
direct evidence is necessary in order to elucidate the true 
origin of the relaxor behavior. For this purpose, neutron 
and X-ray techniques are most suitable, as they can pro- 
vide evidence of local structural ordering. Obtaining 
such evidence is crucial if the formation of polar regions 
is indeed responsible for the relaxor behavior. 

PZN is a prototype relaxor ferroelectric. It was found 
to exhibit a large dielectric dispersion and a broad di- 
electric maximum that depended on both frequency and 
temperature. Earlier studies also reported a structural 
phase transition from cubic to rhombohedral symmetry 
near 413K HU, which falls in the temperature region of 
the maximum of the dielectric peak. More recent stud- 
ies only mention the coexistence of cubic and rhombohe- 
dral phases, with polar nanodomains growing into polar 
microdomains R] and their volume fraction increasing 
with decreasing temperature. Strain appears to play an 
important role in a nanodomain-to-microdomain phase 
transition, which therefore resembles a martensitic phase 



transformation 0. In the last few years, PZN has also 
been intensively investigated in the form of a solid solu- 
tion with PbTi03 ISHlCll , because of its enhanced dielec- 
tric and piezoelectric properties. PZN is already a disor- 
dered perovskite in the cubic phase fll]j , with individual 
Pb atoms shifted from the A-site position along one of 
eight equivalent (111) directions and O atoms forming a 
ring-like distribution within {100} planes. In contrast, 
the Zn/Nb atoms remain located at the normal B-site 
position. In the rhombohedral phase, the displacements 
of Pb and O atoms become correlated, both along a spe- 
cific [111] direction, with Srp^ = 0.31 A and Sro = 0.17 
A at 295K. More specifically, a recent X-ray study |l^] 
of the related system PMN suggests that the three oxy- 
gen atoms may be displaced along three converging [110] 
directions, all three pointing towards the same A posi- 
tion, such that the average oxygen displacement or the 
net oxygen dipole moment in the unit cell is indeed along 
the same [111] direction as that of Pb. It is important 
to note that the displacement of Pb [jllj as well as its 
neutron scattering length are almost twice those of the 
O atom. 

In the present paper, we report the results of a dif- 
fuse neutron scattering study of PZN, in order to iden- 
tify and directly monitor the formation of polar nanore- 
gions resulting from the correlated displacements of the 
Pb and O atoms. This study follows a similar study of 
another mixed perovskite relaxor, Ki-^Li^TaOa (KLT), 
which provided the correlation length or size of the po- 
lar nanoregions as a function of temperature In 
contrast to KLT, PZN exhibits not only substitutional 
chemical disorder but also charge disorder. Moreover, 
due to the lone electron pair effects, Pb 2+ ions form par- 
tially covalent bonds with the neighboring oxygen atoms, 
which results in the off-centering of lead and in an ad- 
ditional local electric dipole moment fbl. Our purpose 



1 



here is to provide direct evidence for these regions and 
to characterize them in a single crystal of PZN. 

In complex perovskites such as PZN, the diffuse scat- 
tering has two possible origins jl5|: i) different scatter- 
ing factors of the two possible atoms (Bi,B2) occupying 
a given B site (j), or ii) static displacements away from 
this site. Both factors arise from concentration inho- 
mogeneities in the compound and may simultaneously 
contribute to the diffuse scattering. It is well known 
that the static structure factor is given by jl6| : 
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where bj is the neutron scattering length, e~ Wj is the 
Debye- Waller factor for the jth atom with coordinate Vj 
in the unit cell and the summation is carried out over all 
positions occupied in the unit cell. The scattering vector 
Q is defined as k — k. In the presence case, small dis- 
placements, 8r, of the atoms from their normal positions 
distort the local cubic symmetry such that rj — (r + 5r)j. 
At high temperature, these displacements vary both in 
magnitude and direction. The scattered intensity is pro- 
portional to F 2 — F.F* so that, experimentally, it is 
possible to determine the amplitude F of the structure 
factor but not its phase. For the sake of argument we 
omit the Debye- Waller factor in the following calculation. 
In accordance with Ref. pH , which reported the Zn/Nb 
atoms at their ideal perovskite position at all temper- 
atures, no displacement is included in cq. (I) for the 
atoms at the B-site. With the origin chosen at the Pb 
site, eq. (1) and for small displacements of the Pb and 
O atoms, the neutron scattering structure factor F can 
then be evaluated at the (1 10) reciprocal lattice point: 

F(Q) = [b Ph + & Zn /Nb - fro ] 

Bragg 

iQ. [Sr Ph b Ph + (Sroi - Sr 2 - Sr 03 ) bo] /rA 
+ v . ' ( 2 ) 

Diffuse 

where 6pb, ^zn/Nb? an d &o are the scattering lengths 
of Pb, Zn/Nb, and O atoms respectively. The scat- 
tering lengths for the Zn/Nb atoms can be found from 
frzn/Nb = g&Zn +§&Nb- In the structure factor above, we 
have separated the Bragg scattering contribution which 
contains information about the average lattice from the 
diffuse scattering which arises from local deviations from 
the ideal structure. These deviations result from cor- 
related displacements of atoms. In the presence case, 
the diffuse scattering is primarily associated with the Pb 
atom, which is strongly off-centered in the PbOi2 unit 
and has almost twice the neutron scattering length of 
the O atom. In addition, because of the Pb displace- 
ment along the [Iff] direction, Jrpjds expected to make 
a greater contribution in that direction than along the 
[001] and [110] directions. 



The neutron scattering experiments were performed on 
the HB-1 triple-axis spectrometer of the High Flux Iso- 
tope Reactor at Oak Ridge National Laboratory. Py- 
rolytic graphite (002) crystals were used for both the 
monochromator and analyzer, to suppress harmonic con- 
tamination. The collimation was 48-20-20-70 and the 
incident neutron energy was 13.6 meV (A=2.453A). Ini- 
tially, the scattering around the (110) reciprocal lattice 
point in the [110]- [001] scattering plane was studied be- 
cause the (110) reflection has the largest structure factor 
0. In each scan, the stepsize was q =0.01 reciprocal 
lattice units (rlu) or 0.016 A^ 1 . The yellowish single 
crystal of Pb(Zni/3Nb 2 /3)03 used in the experiment was 
grown by a high temperature solution growth technique 
fllTf . The measurements were performed upon cooling 
from 550 to 295K and no external electric field was ap- 
plied. 
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FIG. 1. Elastic diffuse scattering of neutrons in the [001] 
and [111] directions near the (110) reciprocal lattice point at 
T=550, 500, 450, 400, and 350K. Inset is scattering intensities 
at T C =413K in the [001]and [111] directions. 

The narrow Bragg peak as well as the relatively broad 
diffuse scattering peak at the (llO)reciprocal lattice point 
are presented in Fig.l for several temperatures. As ex- 
pected, the Bragg scattering is equally strong in both 
the [001] and [111] scans(inset of Fig.l). However, the 
diffuse scattering is more extended in the [001] than in 
the [111] direction, having the form of a rod along a 
cubic direction. It is important to note that the dif- 
fuse scattering is more extended in a direction in which 
the correlation length £ is shorter. We will discuss this 
point later. In order to separate the Bragg and diffuse 
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scattering contributions, the intensity, I(q), was fitted as 
the sum of a Gaussian and a Lorentzian, using a least- 
squares method(Fig.2). In the inset of Fig. 2, the Bragg 
intensity is seen to increase slowly, starting at approxi- 
mately 500K, then more rapidly below 400K. The latter 
rapid increase of the Bragg intensity can be explained 
by i) a relief of extinction due to an increase in mosaic- 
ity below the transition or ii) an unusual behavior of the 
Debye- Waller factor due to freezing of certain atomic mo- 
tion at lower temperatures. A broadening of the Bragg 
peak is normally observed to accompany the freezing pro- 
cess but is not observed here. This would rule out the 
effect of the Debye- Waller factor in the present case and 
favor the first explanation. In fact, a similar behav- 
ior was observed in a single crystal of PMN by Gosula 
et al. Ul using X-ray diffuse scattering. They specu- 
lated that the sharp increase in intensity arose from the 
anti-ferrodistortive ordering that developed at low tem- 
peratures. This anti-ferrodistortive ordering grows to a 
few hundred angstroms with decreasing temperature but 
never reaches truly long-range order. 
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The Lorentzian function used in the fit of the experimen- 
tal I(q) curves, yielded a peak intensity and a full width 
at half- maximum (Aq FWHM ). Such a Lorentzian line- 
shape is predicted by the Ornstein-Zernike model Jl9| : 
/ ~ 2+ g-2 where q = I is the momentum transfer rel- 
ative to the Q=110 Bragg reflection and £ is the corre- 
lation length. The Ornstein-Zernike model assumes a 

correlation function of the form £— — for the polarization. 
Therefore the diffuse scattering width at half maximum, 
Aq FWHM , provides an estimate of £ or, equivalently, of 
the size of the polar nanoregions:£ = -r — where a 

Qfwhm 

is the lattice parameter. The Lorentzian fit parameters, 

1 (?Lax> A< llwHM> and are sh0Wn in Fi S- 3 ( a ) aild ( b ) 

as a function of temperature. In both the [001], and 
[111] directions, the results indicate that: 
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FIG. 3. Temperature dependence of (a). the Lorentzian 
Peak Intensity and square of the FWHM,A^. WHM , (b).the 
correlation length, £ in the [001] and [111] directions. 



FIG. 2. Elastic neutron scattering in the [001] direction at 
T=550 and 350K. The solid curve shows the scattering profile 
fitted by the sum of a Gaussian and a Lorentzian and the two 
components (dash and dot curves) are shown separately. Inset 
is temperature dependence of Bragg peak intensities in the 
[001] and [111] directions. 

Turning to the diffuse scattering component, in Fig. 2, 
the q or Z-dependence of the scattered intensity is pre- 
sented on a semilogarithmic plot at 550K and 350K. At 
550 K the diffuse scattering intensity is slightly asymmet- 
ric which could be related to rather large deformations 
of the crystal lattice giving rise to Huang scattering [^5| . 



i) the diffuse scattering peak intensity, / (<?) max , in- 
creases as the temperature decreases. It is interesting 
however to note that the peak intensity, as measured in 
the [111] direction, goes through a plateau around the 
transition temperature (T c = 413K) before increasing 
again below 400K. This feature is much less visible in 
the [001] direction, for which the intensity increases con- 
tinuously. Fig. 3 (a) also shows that the strongest diffuse 
scattering is observed in the [111] direction, an observa- 
tion that is consistent with the displacement of the Pb 
atoms producing the local polarization in the rhombo- 
hedral phase. The diffuse scattering intensity is in fact 
generally proportional to (P^ oca i) [p0|] . Consequently, an 
increase of the diffuse scattering intensity with decreas- 
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ing temperature reflects the growth of the polar clusters 
or the net correlated displacements of lead and oxygen 
atoms in a [111] direction as reported in Ref. [jll| for 
PZN. It is also in agreement with the structure factor 
calculation F developed above, which predicts a greater 
magnitude of F.F* along the [111] direction than along 
the others. 

ii) the width squared, &-<f FWHM , initially decreases al- 
most linearly with temperature, but then deviates from a 
straight line at approximately 450K, or 40 degrees above 
the transition (T c = 413K) and subsequently levels off 
at low temperature. The initial linear dependence is the 
manifestation of a Curie- Weiss law. This should be ex- 
pected in a temperature range in which the polarization 
fluctuations are dynamic, since the correlation length 
squared, £ 2 , is proportional to the dielectric constant, 
e, p5] which itself is expected to follow a Curie- Weiss 
law initially: 

1 _ £ = c 

A llwHM T-T c 

where C is the Curie constant. The departure of 
Ag^ WHM from a linear dependence, seen in Fig. 3 (a), 
marks the appearance of long-lived polar fluctuations or 
permanent polar regions in the crystal, accompanied by 
local strain fields. This also explains the simultaneous 
increase in Bragg intensity noted above, due to an in- 
crease in mosaicity. It is worth noting that a very sim- 
ilar evolution has been reported by Iwase et al. [jll| for 
the lattice constant, but with a departure from a lin- 
ear dependence at approximately 550K instead of 450K. 
However, their measurements were made on a powder 
obtained from a ground crystal and it is well known that 
relaxors are very susceptible to internal strains. The 
correlation length, £, shown in Fig. 3 (b) increases contin- 
uously with decreasing temperature and saturates at low 
temperatures. This is uncommon for conventional ferro- 
electrics, in which the correlation length decreases below 
the transition. The present behavior indicates that the 
low temperature state is not a homogeneous ferroelectric 
state but, rather, a mesoscopic one, presumably made of 
large polar regions. At T c — 413K, the finite size of the 
polar regions is ~79A (19 unit cells) and ~114A (28 unit 
cells) respectively in the [001] and [111] directions. The 
longer correlation length or effective size of the polar re- 
gions in the [111] direction is found to be £[m] ~ l-2£[ooi] 
(see dotted line in Fig. 3b), consistent with the larger Pb 
displacement component in the [111] direction. This 
result suggests that the polar regions are in the form of 
ellipsoids, preferentially oriented in [111] directions. An- 
other noticeable difference between the two directions is 
the appearance of two kinks in the £[m] curve, at ~ 450K 
and ~ 413K, that are not visible in the £[ooi] curve. As 
discussed above, the first one, at 450K, reveals the rela- 
tively sudden formation of permanent polar regions, and 
the second one corresponds to the transition at T c . The 



absence of these kinks in the £[ooi] curve indicates that 
the growth of the correlation length in that direction is 
only indirectly affected by the formation of the polar re- 
gions at 450K and by their collective realignment at the 
transition at 413K. This result is also consistent with 
the primary net displacements of Pb and three O atoms 
in one of eight possible [111] directions. 

In conclusion, thermal diffuse neutron scattering mea- 
surements in the paraelectric and relaxor temperature 
regions of the PZN crystal clearly reveal the relatively 
abrupt formation of polar nanoregions at 450K, approx- 
imately 40 degrees above the transition at T c — 413K. 
With decreasing temperature, these polar regions grow 
and so does their volume fraction. The correlation 
lengths in the two directions studied are found to be in 
the ratio 6m] ~ 1-2£[qoi] over the whole temperature 
range, except for two kinks in the £[111] curve at 450K 
and 413K. The fact that the correlation lengths remain 
high below the transition indicates that, at low temper- 
ature, the polar regions subsist at low temperature in 
what can be called a mesoscopic ferroelectric phase. Dif- 
fuse scattering results at other reciprocal lattice points 
are presently being analyzed and will be reported in a 
subsequent paper. 
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